Observation of X-ray lines 
from a Gamma- Ray Burst (GRB991216): 
Evidence of Moving Ejecta from the Progenitor 
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We report on the discovery of two emission features observed in the X-ray 
spectrum of the afterglow of the gamma-ray burst (GRB) of 16 Dec. 1999 by 
the Chandra X-Ray Observatory. These features are identified with the Ly Q line 
and the narrow recombination continuum by hydrogenic ions of iron at a redshift 
z = 1.00±0.02, providing an unambiguous measurement of the distance of a GRB. 
Line width and intensity imply that the progenitor of the GRB was a massive 
star system that ejected, before the GRB event, « O.O1M of iron at a velocity 
~ 0.1c, probably by a supernova explosion. 

The nature of the progenitors of GRB's is an unsettled issue of extreme importance 0. 
The merging of a binary system of compact objects (such as black holes, neutron stars, and 
white dwarfs) or the collapse of a massive star (hypernova or collapsar) could all deliver the 
energy required by a GRB, but observational evidence discriminating against the various 
models is still missing. This evidence can be gathered through the measurement of lines 
produced by the medium surrounding the GRB ^|, |3], ffl. However, while observations of 
GRB afterglows have provided much information on the broad band spectral continuum 
and its origin |)[ ||, they have not yet given results of comparable importance on spectral 
lines. In the optical range current measurements are inconclusive, because all of the spectral 
emission lines observed so far are produced by the host galaxy rather than at the burst 
site. The x-ray range appears more promising because theoretical computations show that 
only a dense, massive medium close to the GRB site - such as that expected in the case 
of a massive progenitor - could produce an iron emission line detectable with current x-ray 
instrumentation [7-10]. Indeed, marginal evidence of iron features has been claimed in two 
x-ray afterglows |TT], |T2f but the case is still controversial not only for the limited statistical 



weight, but also for the tight upper limits measured in other afterglows |L3|, and for the 
claimed inconsistency between the redshift derived for GB970828 in x-rays and from the 
host galaxy fI3 ], that could be reconciled only assuming different physical conditions in 
the two bursts |1. The prospect of gathering unique data on the nature of the progenitor 
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made the search for spectral features one of the primary objectives of a Chandra [15[] GRB 
observation program. 

The first Chandra observation of a GRB was performed on the event of 16 Dec. 1999, one 
of the brightest GRB ever detected by the Burst And Transient Source Experiment (BATSE) 
on board of the Compton Gamma-Ray Observatory, with fluence 5 7 > 2.5 x 10~ 4 erg cm -2 
above 20 keV |yj. Following the localization of a strong x-ray afterglow by Rossi X-ray 
Transient Explorer (guided by the rapid BATSE GRB localization) and the characterization 
of its temporal behaviour |fL7ll , and a confirming localization by the interplanetary network 



18] , we estimated that the X-ray flux would only decay to ~ 10 erg cm s by the time 



Chandra could be re-oriented to point at it. This flux level is high enough to employ the 
gratings in conjunction with the Advanced CCD Imaging Spectrometer in the Spectroscopic 
configuration (ACIS-S), and we selected this instrument configuration for the observation. 
Chandra acquired the target on 18 Dec, 04:38 U.T., i.e. 37 hours after the GRB, and 
observed it for 3.4 hours. We found a bright x-ray source |19| with a position (right ascension 
(2000) = 05 h 09 m 31 s .35, declination (2000) = 11°17'05".7) coincident within the 1.5" error 
with the optical 0] and radio |21| transients and with a flux consistent with that expected 
from the XTE extrapolation. 

The spectrum of the x-ray afterglow (Fig. 1) shows an emission line at energy E = 
3.49 ± 0.06 keV. Due to the ubiquity and prominence of iron lines in astrophysical objects 



22] we argue that this line is associated with emission from iron. Some ambiguity remains 



in the rest energy of the emission. Iron K a lines have rest energies ranging from 6.4 keV 
(fluorescence of neutral atoms) to 6.7 keV (He-like ions) or 6.97 keV (H-like ions). In those 
three cases we would obtain redshifts of z = 0.83 ±0.02, z = 0.92 ±0.02 and z = 1.00 ±0.02, 
respectively. In particular, at higher energies the ACIS-S spectrum shows evidence of a 
recombination edge in emission at E = 4.4 ± 0.5 keV. Identifying this feature with the iron 
recombination edge with rest energy of 9.28 keV gives z = 1.11 ± 0.11, consistent with the 
highest of the redshifts implied by the emission line. 
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An iron recombination edge at 9.28 keV is indeed expected when the iron emission is 
driven by photoionization and the medium is heavily ionized by the radiation produced 
by the GRB and its afterglow [fj, |9], |23j. If the medium lies in the line of sight, the edge is 
expected to be seen in absorption at early times |3], |§, and evidence of such a feature has been 
found in another GRB pl[ . At later times, when the medium becomes heavily ionized and 



recombination takes place, the edge is seen in emission. This is our case. In this condition 
x-ray lines are produced almost exclusively through recombination of electrons on H-like 
iron ||25|1 . The measured intensities of the two features are also consistent with theoretical 



expectations (Iedge/^L ~ 0.93(kT/keV) - 2 [|3|], where I e d ge and II are the intensities of the 
recombination edge and emission lines, respectively, and T is the electron temperature of the 
gas). We therefore conclude that the redshift of the GRB is z = 1.00 ± 0.02. We stress that 
this measurement is consistent with the most distant absorption system (z=1.02) found in 



the line of sight towards GRB991216 by optical spectroscopy [p6 |. This system should then 
be in the host galaxy of the GRB, which has probably been identified in deep optical images 

The detection of the line, the measurement of the distance (D = 4.7 Gpc, assuming 

H = 75km s _1 Mpc _1 and q = 0.5) and the fact that the driving process is recombination 

allow us to derive a lower limit on the mass of the line emitting medium. The number of 

iron atoms N Fe needed to produce the observed photon line luminosity^ L = 10 52 L^2 = 

8 x 10 52 photons s _1 is N Fe = L L t rec , where each of the N Fe iron atom produces t~} c 

line photons per second. The recombination time of iron |Q] is t rec — 30Ty s, where 

T = 10 7 T 7 K and n = 10 10 niocm -3 are respectively the temperature and density of the 

electrons. The temperature is constrained from the width of the recombination edge to be 

kT > 1 keV, therefore implying t rec > SOn^s. The total mass of material in the line region 

can be written 

M = M Fe /{X Fe 1.8 x 10- 3 ) > 7X F lL b2 n^M Q (1). 
1 hereafter some quantities are expressed as Y = Y n x 10" 
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Xp e is the iron abundance relative to the sun, where the iron is a fraction 1.8 x 10~ 3 of 
the total mass. The requirement that the ionization parameter £ be high enough to keep the 
Fe in a H-like state, i.e. £ = 4irD 2 F x /nR 2 > 10 4 allows us to estimate nio < 15 where 



we have used the flux observed by Chandra of Fx = 2.3 x 10~ 12 erg cm' 2 s _1 and assumed 
that the distance of the gas from the GRB is R > 2 x 10 15 cm. The latter condition derives 
from the presence of the line 1.5 days after the GRB which, in the cosmo logical rest frame 
of the burst, is (1 + z)~ l times shorter than observed. Using this limit on the density and 
the observed line luminosity allows us to set a lower limit on the mass of: 



M > 5XplM G ; M Fe > 0.01M, 







This large mass is not ejected during the GRB explosion, but in an earlier phase. This 
is the only possible condition to have the material moving at subrelativistic speed (as shown 
below) and be illuminated by GRB photons. The large mass of pre-ejected material excludes 
progenitor models based on double neutron stars, black hole - neutron star and black hole - 
white dwarf. These systems eject material long before they actually merge and the progen- 
itors of these GRB travel far from their formation sites (and their ejecta) before producing 
GRB 0. Conversely, massive progenitors - which evolve more rapidly - lead to a GRB in a 
mass-rich environment. 

Additional information on the origin of the ejecta is derived from the line width [a^ = 
0.23 ±0.07 keV, Fig.l). Thermal broadening (AE/E = (kT/Am p c 2 ) 1/2 ) is negligible here, 
while Compton broadening would require a Thompson optical depth r > 1, which is ex- 
cluded by 0. The line width is therefore kinematical, with v ~ 0.1c. Normal winds 
from stars are not compatible with the parameters of the medium. In fact, the wind 
density is n « 10 3 (M_4i?f 6 2 (w/0.1c) _1 cm -3 . Even for a high value of the mass loss rate 
M_4 = M/(lO _4 M0yr _1 ) = 1, the density is orders of magnitudes lower than that required 
to produce the observed line flux. Weth et al B argued that high density clouds could 
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be produced from a low-density and low-velocity wind, but it remains to be assessed how 
those clouds could be accelerated to high velocities. Alternatively, nearly all the formation 
scenarios of GRB progenitors involve substantial mass loss when the system is in a common- 
envelope phase, a process which is likely to form a diskQ. Interaction of the expanding shell 
of the GRB with the disk could produce a shock-heated gas with density and velocity of the 
same order of magnitude as needed here [10|. The emission from this region can be repre- 
sented by pure thermal plasma, i.e. in thermal and collisional ionization equilibrium, but in 
this case emission from a recombination edge is negligible In fact, to effectively ionize 
iron atoms, electrons should have a temperature comparable to the edge energy, therefore 
producing a feature too smeared to be detected. This may be circumvented if the electron 
population of the shock-heated plasma does not reach complete equilibrium, but the con- 
ditions under which this happens have yet to be studied. The simplest explanation of our 
results is a mass ejection by the progenitor with the same velocity implied by the observed 
line width. The ejection should have then occurred ~ R/v = (i.e., a few months) before the 
GRB. 

The distribution of ejecta and the GRB emission are not highly anisotropic. Let AQ be 
the solid angle of the medium illuminated by the GRB and AR its size. The mass contained 
in the line emitting volume is M = AQR 2 ARnm p . Substituting the limit of the mass derived 
in the left hand side of eq(l) and the limits r = ARnaT < 1 and £ > 10 we derive 

Afi/4vr > mXp l e > 0.1 

where the lower limit corresponds to the extreme case of ejecta of pure iron. The GRB 
emission and the distribution of the medium around it cannot therefore deviate substantially 
from isotropy. The lower limit we have derived on the beaming factor [] is marginally com- 
patible with the estimate of [p0[|. We note, however, that our results refer to the emission 



2 The commom-envelope phase happens when the hydrogen envelop of the secondary star expands en- 
gulfing the compact primary (neutron star or black hole). Friction and tidal forces cause the compact object 
to spiral in the giant's core, ejecting the hydrogen envelope preferentially along the orbital plane, forming a 
disk H 

3 The beaming factor AH./4/ir is the fraction of sky over which GRB photons are emitted. It is equal to 
1 if the emission is isotropical and s=s # 2 /4 in the case of a jet with opening angle 6 
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inside the line emitting region while the measurements mentioned above are based on the 
appearence of a break in the light curve 2-5 days after the GRB, i.e., when the fireball has 
overcome the line region. It is then possible that the initial GRB emission is isotropic and 
it is then collimated by the interaction with a funnel-like medium. This limit on anisotropy 
allows the determination of the total, isotropic, electromagnetic energy produced by a GRB, 
which for GRB991216 is E > 7.2 x 10 52 erg s -1 = 0.04M Q c 2 . Another important implication 
of the previous limit is on the iron abundance of the medium, that has to be much higher 
than solar (Xp e > 60). This high value of the iron abundance indicates that the ejecta were 



- at some stage of the progenitor evolution - produced by a supernova explosion fl29| , |30 . 

In conclusion, the most straightforward scenario that emerges from all the pieces of 
evidence we have gathered is the following. A massive progenitor - like a hypernova or a 



collapsar fl3~l] , |32|j - ejects, shortly before the GRB, a substantial fraction of its mass. This 



event is similar to a supernova explosion, like in the case of the SupraNova model p3 |. 
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Figure 1 The X-ray afterglow spectrum of GRB991216 obtained with the Chandra 
high-energy gratings (High Energy (HE) and Medium Energy (ME) summed together. The 
background is negligible. The exposure time of the observation was 9700 s. In order to 
increase the statistics, the grating spectrum has a bin size of 0.25 A, including about 10 
resolution elements of the ME and 20 of the HE. The dashed line represents the best fit power 
law on the the 0th order ACIS-S spectrum. The peak (i.e. 2 bins) around 3.5 A(E=3.5 keV) 
is detected with 4.7 o confidence. We have verified the robustness of the detection against 
the continuum level. The significance remains above 4<r even when assuming a worst case 



systematic uncertainty in the cross-calibration of the two instruments of 15% [34|. In the 
inset the region on the line is shown with a finer binning. The dotted line represents the best 
fit continuum model to the 0th order ACIS-S spectrum after the addition of a recombination 
edge in emission (see Fig. 2). The line parameters (errors on best fit parameters correspond 
to 90% confidence level for 1 parameter of interest) are II = (3.2 ± 0.8) 10~ 5 cm~ 2 s _1 
E.W. = (0.5±.013) keV, width(a L ) = (0.23±0.07) keV, E = (3.49±0.06) keV. The spectrum 
has been examined at higher resolution to confirm the line broadening. Since each of the 
spectral bins in the figure includes several resolution elements of the instrument, a narrow 
feature would appear in no more than 1 single bin, regardless of how fine is the binning, 
while this is not the case. Deviations around 7 A are ~ 3a, and it is worth noticing that 
they are close to the expected energy (at z=1.0) of the recombination edge of hydrogen-like 
Sulphur. Deviations at ~ 4.4 A are less than 3 o. 

Figure 2 The X-ray afterglow spectrum of GRB991216 obtained with the Chandra 
ACIS-S (0th order). The energy resolution of ACIS-S is 0.1 keV (FWHM) at 4 keV and the 
background is negligible on the whole energy range. The better response at high energies 
compared to the gratings allows us to single out the presence of a further emission feature. 
Fitting a model (continuous green line) composed by an emission edge (blue dashed line) plus 
a power law (green dashed line) plus line ( orange dashed line) provides a satisfactory fit to 
the data (xl — 0.95, v = 26). The addition of the edge improves the fit by Ax 2 /xl — 16.3, 
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that corresponds to a confidence level of 99.5% (F test). Best fit parameters of the edge are 
E = AA ± 0.5 keV, I edge = (3.8 ± 2.0) lO^cm" 2 s^ 1 and width a edge > 1 keV. For the power 
law we derive (F x (2 - lOkeV) = 2.3 10~ 12 erg cm" 2 s _1 , photon index T = 2.2 ± 0.2, N H = 
(0.35±0.15) 10 22 cm -2 , consistent with the absorption in our Galaxy (N HG = 0.21 10 22 cm~ 2 ). 
Line parameters are consistent with those derived from the grating. Moreover, the edge is 
consistent with the grating data, as shown by the dotted line in the inset of Fig.l. 
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